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We show that the decorated honeycomb ( "star" ) lattice supports a number of topological insu- 
lating phases with a non-trivial Z2 invariant and time-reversal symmetry protected gapless edge 
modes. We investigate the stability of these phases with respect to various symmetry breaking per- 
turbations and demonstrate the connection to the recent discovery of an exactly solvable S — 1/2 
chiral spin liquid model [Phys. Rev. Lett. 99, 247203 (2007)] with non-Abelian and Abelian ex- 
citations on the same lattice at strong interaction strength. Our work highlights the relationship 
between non-interacting topological band insulators and strongly interacting topologically ordered 
spin systems, and points to promising avenues for enlarging the number of known examples of both. 
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I. INTRODUCTION 

In recent years, the study of various types of topo- 
logical order in condensed matter physics has dramat- 
ically increased. 1-3 The interest in this topic has been 
driven in large part by the fractional quantum Hall ef- 
fect, and efforts to understand the high temperature su- 
perconductors. In both cases electron interactions are 
fundamental to the phenomena. However, a new class of 
systems, non-interacting Z2 topological band insulators 
(TBI) with time-reversal symmetry (TRS), has diverted 
attention to topological properties that do not depend on 
interactions (but are robust to weak interactions). 4 ' 5 The 
existence of topological properties in models that can be 
exactly solved in the non-interacting limit, and treated 
to a high degree of accuracy by conventional band the- 
ory methods in the interacting limit, has led to precise 
predictions for experiment. 6-9 An unusually rapid verifi- 
cation of many of these predictions in experiment has fol- 
lowed, and there are now several known examples of this 
state of matter in both two dimensional 10,11 and three 
dimensional systems 12-14 . In some of these materials, 
topological properties are expected to be robust up to 
room temperature and therefore hold great promise as 
components of future electronic devices. 12 

Current theoretical research on topological insulators 
is proceeding along several parallel tracks. On the one 
hand, there is great interest in identifying new phys- 
ical systems that will possess topologically non-trivial 
phases 15-21 , while on the other hand there are fundamen- 
tal questions about the fate of topological properties as 
the strong electron interaction limit is approached 22 ' 23 . 
In this work, we contribute to both directions by provid- 
ing several other examples of Z 2 TBIs on a lattice where 
they have not been reported before-the decorated honey- 
comb lattice. We also establish a topological connection 
at 1/2 filling between the non-interacting limit and the 
strongly interacting limit where an exactly solvable elec- 
tron model (the Kitaev spin model) is realized on the 
same lattice. 24 We are unaware of any other model that 
realizes exactly solvable states at weak and strong inter- 
action, both with topological properties. Moreover, via 



explicit calculation, we show these two limits share topo- 
logical properties, even though their symmetries are very 
different. 

Our discussion focuses on a tight-binding model of 
fermions hopping on the 2-dimensional decorated hon- 
eycomb lattice shown in Fig. 1. This lattice is a "cousin" 
of both the honeycomb lattice and the kagome lattice, 
each of which is known to support TBI phases. 4 ' 15 ' 25-27 
In a certain regard, the decorated honeycomb lattice can 
be viewed as an "interpolating" lattice between the hon- 
eycomb and the kagome: If one shrinks the triangles at 
the verticies of the underlying honeycomb lattice (sites 
with hopping parameter t in the figure) to their cen- 
ter point, the honeycomb lattice is recovered, while ex- 
panding the triangles until their corners touch produces 
the kagome lattice. One might consider this geometrical 
property to be the key reason the decorated honeycomb 
lattice supports topological insulator phases, given that 
the honeycomb and kagome lattices also support topo- 




FIG. 1: (color online) (a) The decorated honeycomb lat- 
tice has a triangle at each vertex of the honeycomb lattice. 
The 6-site unit cell with "sublattice" A and B is contained in 
the parallelogram indicated by ai and a2. Nearest neighbor 
hopping on vertex triangles occurs with amplitude t, between 
triangles with amplitude t' , and with ±iAso for second neigh- 
bor hopping as indicated in (b). Topological phases occur at 
a number of filling fractions (see Fig. 3), as well as in the case 
that Aso = and t is allowed to be complex corresponding 
to finite flux through vertex triangles, (c) The first Brillouin 
zone including a path along the high-symmetry lines. 
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logical phases. However, because the unit cell of the dec- 
orated honeycomb lattice contains 6 sites (compared to 
3 for the kagome and 2 for the honeycomb) its phase di- 
agram is much richer than that of either of its "cousins" 
and some novel features appear that we will discuss in 
more detail below. 



II. TIGHT-BINDING MODEL 

The Hamiltonian for our (initially non-interacting) 
problem is 

H = H a + H so + i?cD\v + H R - (1) 
The nearest-neighbor hopping is described by 

H = -tJ2 4<> - *' S 4«7a + h.c. (2) 
(ij),<r,A (y),<7,A->A 

with amplitude t on the triangles "A" and with ampli- 
tude if between triangles "A — > A" , as shown in Fig. 1 (a) . 
The intrinsic spin-orbit coupling, 

Hso = «A so ^ eij ■ s a pcl a Cj0 + h.c, (3) 

describes the second-neighbor hopping with amplitude 
±i\so, see Fig. 1(b). The sign of the amplitude is differ- 
ent for different spin orientations s z — ±1, sis the vector 
of Pauli matricies and = (d^ x d 2 .)/|d^ x d?-| is a 
vector normal to the x — y plane describing how the path 
((ij)) was traversed using the standard conventions. 4 In 
Eq. (1), -ffcDW an d are charge density wave and 
Rashba spin-orbit terms, respectively. The CDW Hamil- 
tonian is 

#CDW = X] ^viCiaCia, (4) 

and the Rashba Hamiltonian is 

Hb. = «Ar ^ c L(^«/3 X d ij) z Cjp + h.c, (5) 

where A„j is an on-site potential possibly differing on each 
of the 6 unit cell sites labeled in Fig. 1, Ar is the strength 
of the Rashba coupling and d%j is the unit vector connect- 
ing site i to j. 

III. PHASE DIAGRAMS 

The 6 (doubly degenerate) bands coming from the 6- 
site unit cell (see Fig. 1) for H Q + H S q are shown in 
Fig. 2 along the various high symmetry directions. The 
first Brillouin zone is identical to that of the honeycomb 
and kagome lattices which share the same underlying tri- 
angular Bravais lattice, 4,15 see Fig. 1(c). There are Dirac 
points at K and K', two quadratic band crossing points 
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FIG. 2: (a) and (b) show the band structure of the tight 
binding model Ho + Hso with t = t' along the path shown in 
Fig. 1(c). In (a) Aso = and in (b) Aso = O.lt. There are 
Dirac points at K and K' (not shown) and quadratic band 
crossing points (QBCP) at V in (a), while in (b) Aso 7^ 
opens up a gap at each of these points and destroys the flat 
bands, (c) The phase diagram at / = 1/2 (involving QBCP) 
with Aso = O.lt. (d) The phase diagram at / = 1/6 (involving 
Dirac points) for Aso = O.lt. We have chosen a staggered 
sublattice potential configuration where all the sites in A(B)- 
triangle (see Fig. 1) have potentials A„ (— X v ). 



(QBCP) at T, and two flat bands present when Aso = 0. 
We note that the lower QBCP appears at filling fraction 
/ = 1/2 for if < it/2 and at / = 1/3 for if > it/2. Sim- 
ilar band features are also found on the kagome lattice 
at the same Brillouin zone points. 15 ' 27 When the second 
neighbor hopping Ago 7^ 0: a gap opens at the Dirac 
and the QBCP and topologically non-trivial phases ap- 
pear; denoted as quantum spin Hall (QSH) insulator in 
Figs. 2 and 3. By explicitly computing the invariant 
using the parity eigenvalues at the time-reversal invari- 
ant moment 6 and checking for helical edge states in a 
strip geometry 4 , we have found the phase diagrams for 
different filling fractions, /. The results are summarized 
in Fig. 3. 

One feature of the decorated honeycomb lattice that 
differs from the kagome and honeycomb lattices is the 
natural presence of two (t and if), rather than one, 
nearest-neighbor hopping parameters. This effectively 
adds an additional degree of freedom to the phase di- 
agram and can lead to transitions to topologically non- 
trivial phases even when there is not an obvious Dirac 
point or quadratic band crossing involved in the nearest- 
neighbor hopping model, such as occur at filling fraction 
1/3 in Fig. 2(a). As Fig. 2(b) shows, when Aso is turned 
on, an "incipient" band touching point develops at the T 
point for filling fraction 1/3 and this effectively drives the 
transition to the topologically non-trivial state. Thus, 
the band structures with zero spin-orbit coupling do not 
always clearly reveal potential topological transitions for 
strong spin-orbit coupling. 



3 




M 



-2 2 

Aso/t 




-2 2 

AsoA 


: f 


2 

SO 

> 




2 2/3 i 


-2 


5/6 




-2 2 

AsoA 




-2 2 

Aso/i 



~T o 

~ 2 |l/2# \ 
-2 2 

Aso/t 

U QSH insulator 
I I band insulator 
I I semimetal 



FIG. 3: Phase diagrams for the decorated honeycomb lattice 
with t and t' real in the absence of a staggered on-site poten- 
tial and no Rashba coupling. Several filling fractions / are 
shown (lower left corner). For fixed / and Aso it is possible 
to drive a transition between a topological insulator and a 
non-topological phase by varying the ratio t' /t. 



At / = 1/3, 1/2, and 2/3 there are (electron-hole com- 
pensated) intervening metallic phases between topologi- 
cally trivial and non-trivial insulators. This also indicates 
that a "direct" transition coming from a band inversion 
is not generic in this model. 28 The filling fractions with 
such an intervening metallic phase mimic the behavior 
of disorder on the honeycomb lattice in the presence of 
finite Rashba coupling. 29 Also note that varying the ratio 
of t' ft at fixed Aso can lead to a transition between a TBI 
and a trivial insulator. More surprisingly, increasing Aso 
for fixed value of t' /t can trigger a transition between a 
TBI and a trivial insulator, as seen for / = 1/3 and 1/2. 

Next we turn to an analysis of the stability of the 
topological phases indicated in Fig. 3 in the presence of 
Rashba interaction and on-site (CDW) potentials. The 
stability of the topological phases at / = 1/6 and / = 1/2 
is shown in Fig. 2 (c) and (d), where we used A„, = A„ 
on "sublattice" A and X V i = —X v on "sublattice" B as 
shown in Fig. 1. Consequently, the stability regions are 
qualitatively similar to the analagous model on the hon- 
eycomb lattice. 4 ' 30 The stability region of the QBCP at 
/ = 1/2 is larger than that for the Dirac point at / = l/6, 
which we attribute to a larger value of the gap (~3 times) 
at the QBCP when \ R = X v = 0. 31 



IV. EFFECTIVE LOW ENERGY DESCRIPTION 
AT DIRAC POINTS 

At the Dirac points it is straightforward to derive an 
effective low energy description for arbitrary X V i in the 
6-site unit cell. Setting the zero of energy to be right at 
the Dirac point (either at / = 1/6 or / = 2/3) for Aso = 
A_r = A„j = 0, the effective low-energy Hamiltonian is 



given by 



H' = H' + iJg 



with 



H 'so 

^CDW 



avFh(k x r z a x + k y <Ty), 
-4aw(t'/t) X so <t z t z s z , 
-a w R (t'/t) X R (a x T z s y - a y s x ), 
g I + a(g x a x + g y T z a y + g z cr z ). 



(6) 



(7) 
(8) 
(9) 
(10) 



We have adopted aT z ,Ci,s, notation similar to Ref. [32]: 
The r z = ±1 describes states at either the K or K' 
points, the <y z = ±1 describes the two bands that are 
involved in the Dirac band crossing (analog of A and B 
sublattice bands on the honeycomb lattice), and s z = ±1 
represents the electron spin as it did in Eq. (3) and 
Eq. (5). The parameter a — ±1 refers to the Dirac point 
at / = 1/6, and / = 2/3, respectively. We have also 
defined two functions describing the dependence of the 
effective low-energy theory on the parameter x = t' /t: 



w(x) 
Wr{x) 



2\/9~+4r 2 ' 
3 + 2V3x 



4x 2 



The effective Fermi velocity entering Eq. (7) is 

vp = w(t' /t)vo 

where v = ta/h and a is the length of the unit cell 
vector. It follows from Eq. (9) that for A OT ; = Ar = 
the spin-orbit coupling opens up a gap with magnitude 
-Egap = 8\w(t' /t)Xso\. The parameters entering the low- 
energy description of the CDW term, Eq. (10), are given 
by 
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9x 

9y 

9z 



X v \ + X V 2 + A„3 + A„4 + A„5 + A„6 

A„i + A„2 — 2 At,3 + A„4 + A„5 - 
3^ 

—X v i + X V 2 — A„4 + A„5 



w{t'/t) 
w(t'/t) 



2A 
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w(t'/t) 



~ Aul — X V 2 — A t ,3 + A„4 + A„5 + A^6 

vWJt) ' 



A finite go can be absorbed in a shift of the chemical 
potential. 

It is useful to consider a few important limits of the 
general low-energy form of i?cDW- First take X V i = X v 
for sites on the A-triangle and X vi = — X v for sites on the 
B-triangle. In this case, 

-Wcdw — fJz ' 

which is identical to the form of the expression for the 
honeycomb lattice and will generically open a gap at the 
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Dirac point. 4 We have verified that the low-energy de- 
scription given above produces the same stability phase 
diagram and phase boundary shown in Fig. 2(d) as a di- 
rect diagonalization of the full 6-band Hamiltonian. An- 
other important limit to consider is that of general A„j. In 
that case, the physics more closely resembles the kagome 
lattice where an effective axial gauge held appears 15 with 



At, i + \ V 2 — 2 At,3 + A„4 + A U 5 — 2 A „g 
3^ 



I, 



• i _ A^i — A„2 + A„4 — At,5 



when A^i + A„2 + A.„3 — A„4 — A„5 — A„e = 0, where / = ±1 
refers to the two Dirac points K and K' . If, on the other 
hand, A^i+A^+A^-A^-A^-A^ ^ a gap 



E, 



gap 



|Aj,i + A„2 + A„3 — ^v4 ~ A„5 — A^6 1 

v/9 + 4(i'A) 2 



opens with a smallest direct gap at shift A l x ,A l y with re- 
spect to K or K' . Thus, the behavior of the decorated 
honeycomb lattice with respect to Hqdw is another ex- 
ample of the ways in which this lattice "interpolates" 
between the honeycomb and kagome lattices, and we ex- 
pect, for example, analogs of the kekule phase to be re- 
alized as well. 15 ' 33 



V. ADIABATIC DEFORMATIONS AND THE 
KITAEV MODEL 

We now turn our attention to one of the features of the 
decorated honeycomb lattice which is related to its geom- 
etry: Topological phases exist even in the absence of sec- 
ond neighbor hopping when t is made complex (obtained 
by putting a flux through the vertex triangles) . Below we 
show by an explicit calculation for spinless fermions that 
the model obtained in the absence of second neighbor 
hopping but with complex t (and possibly also complex 
if) can be adiabatically deformed into a model with real 
t, if and Ago- An example of such an adiabatic defor- 
mation is illustrated in Fig. 4 and in the last part of this 
section we will describe each step of the deformation in 
detail. 

The adiabatic connection we establish also holds for 
time reversal invariant models of electrons with spin: For 
s z conserving models on the honeycomb lattice, Kane and 
Mele showed 4,32 that one can view a Z 2 TBI in 2-d as 
two copies of Haldane's model 26 with different effective 
magnetic fluxes (with a net zero flux through the unit 
cell) for different spins (so that under time-reversal each 
copy is transformed into the other, thus preserving TRS 
overall). Moreover, as long as the gap does not close, 
s z non-conserving terms are also allowed. With this in- 
sight, it is evident that any lattice model that supports a 
quantum Hall effect for spinless fermions will support a 
Z 2 TBI for electrons with spin (by taking the appropriate 
"second copy"). 
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FIG. 4: (color online) Schematic illustration of the continu- 
ous path which adiabatically connects the model (o) and (i) 
lacking second neighbor hopping but having complex t (the 
representative free fermion model of the ground state sector of 
the Kitaev model on the decorated Honeycomb lattice) with 
the spinless model (v) with real t, if and Aso, denoted by 
(v). Along this path, the continuous deformation does not 
lead to a gap closing and the Chern number stays constant. 
This establishes the topological connection between the two 
models. 



There is an interesting consequence of the above men- 
tioned adiabatic connection. It allows us to topologically 
relate the phases in the non-interacting tight-binding 
model at half-filling, see Fig. 3, to the chiral spin liq- 
uid phases recently reported in the Kitaev model 24,34 on 
the same lattice which can be viewed as a strongly in- 
teracting electron model with spin-orbit coupling 35 . The 
Kitaev model on the decorated honeycomb lattice is de- 
fined in the following way 24,36 



H = 



;r-link 



• E VM 

y-link 



z-link 



+ e + E ■'»". + E J '*°w 

y'-link 



:r'-link 



z'-link 



In the summation, i and j are nearest neighboring sites 
connected by a a-link (a = x,y, z,x' ,y' , z') as shown 
in Fig. 1(a). After a Jordan- Wigner transformation 
this model can be mapped onto a model of free majo- 
rana fermions hopping in the background of static Z 2 
fluxes. The ground state spontaneously breaks time- 
reversal symmetry and is described by a chiral spin liquid 
with either Abelian or non-Abelian vortex excitations. 24 
In the following we set J x = J y = J z = J > and 
J' x = J' y = J' z = J' > 0. The ground-state sector cor- 
responds to a uniform flux configuration. It possesses a 
representative free fermion model 25 



ffcsL = 7E<X°>)^ 



kfj 



(11) 



and 7 = ±1 specifies the way the time reversal symmetry 
is spontaneously broken. Here, ip^l' a = l>--->6 are 
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fermionic annihilation (creation) operators and we have 
defined the matrix 
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The flux pattern derived from the matrix ffw is illus- 
trated in (o) of Fig. 4. The spectrum of Eq. (11) is gapped 
at half filling as long as J' ^ y/3J and the Chern number 
is v = -7 for J' < V3J and v = for J' > V3J. 24 The 
two sectors connected by time reversal symmetry and 
characterized by the parameter 7 = ±1 are similar to 
the two "copies" characterized by s z = ±1 in the Kane- 
Mele type model. However, in the strongly interacting 
limit (Kitaev model) the system spontaneously chooses 
one sector whereas the non-interacting TBI model in- 
volves a summation over the two sectors (spin). 

We now discuss the adiabatic connection of the model 
H^°' to the spinlcss model with real t, t' and Aso; 
model in Fig. 4. As the starting point for the 

continuous deformation we use a gauge equivalent pat- 
tern (i) obtained from (o) by replacing -0k4 — > —^4 and 
V J k5 ~ "0k5- From (i) to (ii) the phase difference for 
hopping between the triangles is gradually reduced to 
zero. This can be achieved by replacing ±i J' in by 
exp(±is7r)J' and continuously reducing s from 1 to 0. 
This process does not change the fluxes through the tri- 
angles and the dodecagons; instead the global fluxes are 
modified and the whole spectrum is moved in k-space 
according to k\ — > fcj + (1 — s)n and k 2 —> k 2 . Conse- 
quently, the direct gap stays constant. From (ii) to (hi) 
a continuous gauge transformation is applied to change 
the phase difference for hopping within the triangles from 
±7r/2 to 7r/6 (this does not modify any flux and the gap 
remains constant). From (hi) to (v) we turn on Aso a s 
shown in (iv) and make t real. Along this path the value 
of the gap varies in general. Explicitly, we can define the 
set of matrices 



r(<Mso) = 

/ e"^J e^J J'e~ ik2 D -E \ 

e i4>j o e -i<t>j _ D j/ e ;fei F 

e -i4>j e i<l>j E -F J' 

J'e lk2 -D* E* e'^J e^J 

D* J'e- lkl -F* e^J e'^J 

\ -E* F* J' e'^J e^J J 



Here we have introduced 
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FIG. 5: (color online) Contour plot of (a) the Cern number 
and (b) the gap of the model defined by F(0, Aso) at half 
filling for J = J'. Also shown is a possible path which adia- 
batically connects the flux patterns (iii), (iv) and (v) defined 
in Fig. 4. 



Clearly, 1^/6,0) = and r(0,A SO ) = # (v) - By 

numerical examination of the gap and of the Chern 
number 37 at half-filling we find that there is a large range 
of parameters which allows one to adiabatically connect 
the model (iii) with the model (v) . This is shown in Fig. 5 
where we plot the Chern number in (a) and the value of 
the gap in (b) obtained for J = J' . 

The model defined by T(4>, Aso) shows a complex 
phase diagram with a variety of topological phases dis- 
tinguished by different values of the Chern number. As 
long as there is a direct gap, the Chern number is well- 
defined and regions with different values are necessarily 
separated by gap closings. However, there are also re- 
gions in parameter space where an indirect gap is closed 
indicating the presence of partially filled bands at half- 
filling, see Fig. 5(b). The phase diagram of r(<^>, Aso) 
also depends on the ratio J'/ J. For J' > ^/ZJ the Chern 
number of ij( m ) is zero and the connection (iii) to (v) 
holds between topologically trivial phases. 

These arguments demonstrate the adiabatic connec- 
tion between the strongly interacting chiral spin liquid 
phases of the Kitaev model and the phases obtained from 
the spinless model at half-filling with real t, t' and Aso- 
In this sense, we also establish a connection to the Z 2 
TBI when two "copies" of the spinless model are taken. 
It can be shown that a generalized spin-3/2 Kitaev model 
on the kagome lattice also supports a chiral spin liquid 
ground state with non-Abelian excitations 38 and similar 
arguments can be made for connections to other topolog- 
ical phases. 39 



VI. SUMMARY 

In summary, we have shown that the decorated hon- 
eycomb lattice supports Z 2 topological phases at various 



(> 



filling fractions, discussed their stability, and described 
the similarities and differences with the topological 
phases on the kagome 15 and honeycomb lattices 4 . The 
limit of weak spin-orbit coupling yields phase diagrams 
which are very similar to the ones observed on the honey- 
comb or kagome lattice. This observation is in agreement 
with the conclusions drawn from the effective low-energy 
theory at the Dirac points. The situation for strong spin- 
orbit coupling can be quite different and leads to novel 
aspects. One surprising observation is that a large Aso 
can trigger a transition from a TBI to a trivial insulator. 

We have also shown that the tight-binding models with 
real t, t' and Aso can be adiabatically connected to mod- 
els without second-neighbor hopping but with complex t 
(and possibly also complex t'). This property was explic- 
itly demonstrated at half filling by a calculation of the 
gap and the Chern number. Moreover, we have argued 
that this adiabatic connection allows us to topologically 
relate the chiral spin liquid phases recently discovered on 
the decorated honeycomb lattice to the phases obtained 
at half-filling in the non-interacting TBI model. Our 
work therefore provides an example of a non-interacting 
and a strongly interacting model defined on the same lat- 
tice which are both exactly solvable and show topologi- 
cally related states. To determine the precise form of the 



spin-orbit interaction, the number, and the size (coupling 
strength) of the terms needed in a generalized extended 
Hubbard model at half filling to interpolate between the 
topological band insulator and the Kitaev model on the 
decorated honeycomb lattice is an interesting open prob- 
lem beyond the scope of this work. However, based on 
our results here and related studies that realize Kitaev 
models in certain low energy limits 35,40 we believe that 
such an interacting microscopic model can be found. 

Finally, we note that an underlying "star" lattice has 
been experimentally reported for a polymeric Iron(III) 
acetate 41 , and some of our results may be relevant for this 
solid state example of a decorated honeycomb lattice. We 
also believe it is possible to realize much of the discussed 
physics (including Kitaev models), in cold atomic gases, 
given that its two cousins, the honeycomb and kagome 
lattices, can be realized in optical lattices. 42 ' 43 
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